Universal Finite Temperature Properties of a Three Dimensional Quantum 
Antiferromagnet in the Vicinity of a Quantum Critical Point 
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We consider a 3-dimensional quantum antiferromagnet which can be driven through a quantum 
critical point (QCP) by varying a tuning parameter g. Starting from the magnetically ordered 
phase, the Neel temperature will decrease to zero as the QCP is approached. From a generic 
quantum field theory, together with numerical results from a specific microscopic Heisenberg spin 
model, we demonstrate the existence of universal behaviour near the QCP. We compare our results 
with available data for TlCuCl^. 
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The subject of continuous Quantum Phase Transitions 
(QPT's) and the behaviour of quantum systems in the 
vicinity of the corresponding quantum critical points is 
a frontier area of research both in theory and in experi- 
ment. P, Q A QPT is a transition at zero temperature, 
in the nature of the ground state, and is due to quan- 
tum fluctuations that can be enhanced or suppressed by 
varying some coupling constant. In real materials QPT's 
can be driven by pressure, by applied magnetic field, or 
by some other parameter. 

In the present work wc consider an 0(3) QPT which 
occurs between a magnetically ordered Neel phase and 
a magnetically disordered 'valence-bond-solid' (VBS) 
phase in a class of SU(2) invariant Heisenberg spin sys- 
tems. This problem has attracted a great deal of atten- 
tion in recent years, mainly in two-dimensional (2D) sys- 
tems. It has been established that the interplay between 
quantum fluctuations and thermal fluctuations at low but 
finite temperatures influences the dynamics in the vicin- 
ity of a QPT in a highly nontrivial way 0, 1^. However, 
in 2D systems there is no finite temperature magnetic 
order, due to the well known Mermin- Wagner theorem. 
One would expect that in 3D systems (3D -I- time) the 
presence of a finite Neel temperature and an extended re- 
gion of magnetic order will affect the interplay between 
quantum and thermal fluctuations, and lead to new fea- 
tures not seen in 2D. An obvious question is the nature 
of the vanishing of the Neel temperature and its scaling 
with the magnetization and with the coupling constant 
as the QPT is approached. To the best of our knowledge 
the generic problem of the finite temperature behaviour 
of 3D systems in the vicinity of an 0(3) QPT has not 
been previously considered. The present work addresses 
this question. 

Specifically, we discuss three aspects of this question. 
The first is to consider a general Landau-Ginzburg field 
theory, which is independent of the details of any mi- 
croscopic model, and hence generic. The predictions of 
this approach are then compared with experimental re- 
sults for the material TlCuCla. Finally we present re- 
sults obtained for a specific microscopic Heisenberg spin 



model, obtained using a variety of series-expansion meth- 
ods. While the numerical precision close to the QPT is 
only moderate, the results are consistent with the field 
theory predictions, and reinforce our conclusion that the 
behaviour is universal. 

To develop a quantum field theoretic description we 
start from the standard effective Lagrangian describing 
an 0(3) QPT, of the form 



£ = i (<^- [0 X B] 
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In the present work we consider zero magnetic field, 
B = Q. The vector field describes the staggered magne- 
tization. The QPT results from the mass term, assumed 
to be of the form = ^2^^ _ ^^-j^ where > is 
a coefficient and g is a coupling parameter (In TlCuCla 
the coupling parameter is an external hydrostatic pres- 
sure). When g > Qc the mass squared is positive and this 
corresponds to the magnetically disordered phase with 
gapped triply degenerate excitations. These are some- 
times called 'triplons' but we will use the term 'magnon' 
in both phases. The zero temperature gap is 



A = m = g — Qc 
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When g < gc the mass squared is negative and this results 
in a nonzero expectation value 



1(^)1 = 
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that describes the spontaneous staggered magnetization 
at zero temperature. This is a magnetically ordered 
phase with a gapped longitudinal mode and two trans- 
verse gapless Goldstone modes. We note that Lp has di- 
mensions of energy~^^^ , and therefore cannot be directly 
compared with the dimensionless staggered magnetiza- 
tion. The zero temperature energy of the magnetically 
ordered ground state is 
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FIG. 1: The phase diagram in the vicinity of a QPT. The 
blue triangles connected by the blue solid line show results 
of series expansion calculations of the spin wave gap in the 
magnetically disordered VBS phase. The Neel temperature 
is shown a.t g < gc where the system is magnetically ordered 
at T < Tjv. The red squares connected by the red dashed 
line show results of series expansion calculations of Tm- The 
magenta solid line shows the field theory prediction for Tjv. 



The generic phase diagram is shown in Fig. [T] The 
specific parameter values shown in the figure correspond 
to the particular model which we consider below. 

The magnetic ordering at g < gc is destroyed at T > 
T/v. To find T/v we calculate the sclfencrgy E shown in 
Fig.[2l The four- leg vertex in Fig. [2] is due to the quartic 




FIG. 2: Magnon selfenergy 

a-term in ([1]). To calculate the single loop selfenergy in 
the magnetically disordered phase it is sufficient to decou- 
ple the quartic interaction, a[(jp]'^ — ^ a{tfP)i^ — >■ 'EijP. 
When performing the decoupling, one has to be care- 
ful about the combinatorial factor that is due to various 
ways of the field couplings. A straightforward calculation 
gives the following selfenergy in the magnetically disor- 
dered phase {g > gc or T > Tn at 5 < gc), 
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where i is any of three Cartesian components of ip, and 
nui = (flkjOki) = l/(e"''/"^ — 1) is the thermal population 
of this component. The quantum fluctuation part of ([5]) 
is ultraviolet divergent. 
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Here A is an ultraviolet cutoff and 6 is the gap in the 
spectrum, for example at T = and g > gc, 6 = A. The 



quadratically divergent part oc A^ of the self energy has 
to be removed by renormalization. In other words, this 
part is absorbed in the value of the critical coupling con- 
stant gc- The logarithmic part depends on both the ultra- 
violet cutoff A and the infrared cutoff 6, and is therefore 
a real physical correction. However, we expect this loga- 
rithmic correction to be small and therefore we disregard 
it, (see also the discussion in Ref. 0). The parameter 
that suppresses the correction is the prefactor I/tt^, and 
in essence it is related to the 3D character of the prob- 
lem. Neither the existing experimental data presented 
below nor results of numerical simulations also presented 
below have a sufficient accuracy to pin the logarithmic 
corrections down. All in all this implies that the entire 
quantum part of the self energy is renormalized out. 
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where the subscript 'R' stands for 'renormalized'. At 
T = T/v the excitation spectrum is gapless, (5 = 0, cjk = 
ck. Hence a calculation of the integral in ^ gives S^j = 

^"^S ■ If the magnon spectrum is anisotropic with three 
different principal velocities then has to be replaced 
by C1C2C3. The magnon gap at the Neel temperature is 
= + Sr = 0, and hence 



zero. 
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Thus, the Neel temperature is directly proportional to 
the zero temperature staggered magnetization ([3]). A 
similar scaling was obtained recently in Monte Carlo sim- 
ulations with various kinds of models 0] . 

Equation ^ can be compared with experimental 
data for TlCuCla . Values of gaps at zero temperature 
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FIG. 3: Zero temperature magnon gaps and Neel temperature 
in TlCuCls versus pressure. All values are given in meV. The 
system is magnetically disordered at p < pc = 1.07kbar, and 
magnetically ordered at p > pc. Symbols show experimental 
data from Ref. Q| and curves show theoretical results. 

versus pressure are plotted in Fig.|31 The critical pressure 
is pc = l.OTkbar. Note that Fig. [3] is 'mirror reflected' 
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compared to Fig. [1] the magnetieally ordered phase is 
at p > pc, therefore, to compare with Eqs.®,©,© 
we choose g = —p. In the ideal situation correspond- 
ing to the action ([T]) one should expect triply degenerate 
gapped excitations in the magnetically disordered phase 
ai p < Pc, as well as one longitudinal gapped mode (A^) 
and two gapless transverse modes in the magnetically or- 
dered phase at p > Pc- In the real compound there is 
a small easy plane anisotropy and due to the anisotropy 
one of the transverse magnons in the magnetically or- 
dered phase is gapped, with A^, = 0.38meV. For the same 
reason the triple degeneracy &t p < pc is lifted. Disre- 
garding the small anisotropy effects and using Eq.Q we 
fit the gap in the magnetically disordered phase. The fit 
is shown in Fig. [3] by the blue dashed line, and results in 

1/2 

the following value of A, A ~ 0.66 meV/kbar ' . Other 
parameters of the effective action ([T]) were determined 
in the analysis of magnon spectra and Bose condensa- 
tion of magnons performed in Ref. 0, ci = 7.09meV, 
C2 = 1.12meV, cg = 0.51meV, a = 21(1 ± 0.2) meV^. 
Substitution to Eq.(l7|) gives the theoretical prediction of 
the Neel temperature plotted in Fig. [3] by the solid ma- 
genta curve with error bars that are mainly due to uncer- 
tainty in the value of a. This curve is very close to the 
experimental points, shown by magenta stars. The exper- 
imental values of the Neel temperature are slightly higher 
compared to the theory, especially close to the QPT 
point. As one would expect the magnetic anisotropy, 
pointed out above, leads to an enhancement of the Neel 
temperature. We believe this fully accounts for the dis- 
crepancy. 

While the above theory is generic, and independent of 
the details of any microscopic model, it is interesting and 
important to consider a specific model and compare re- 
sults with general theory. A specific microscopic model 
can be analysed only numerically, so below we consider 
a sort of numerical experiment versus the real experi- 
ment discussed above. Many previous numerical studies 
of QPT's have been reported. These have been largely 
based on Heisenberg spin models in which the system 
can be tuned through a QPT by varying a particular 
coupling parameter in the Hamiltonian. Most of these 
models have been two-dimensional. Examples include 
antiferromagnets with strong and weak bonds, with or 
without frustration , and bilayer systems ]ll, [3], 

where the QPT separates a conventional Neel antiferro- 
magnetic phase from a spin-dimerized phase with only 
short range correlations and no magnetic order. 

Here we consider a 3D spin-1/2 Heisenberg antiferro- 
magnet. Our model, shown in Fig|3l^a), has weak and 
strong bonds of strength J and gJ respectively. For 
(7 = 1 we have an isotropic cubic antiferromagnet, which 
has reduced st agg ered magnetization in the ground state 
(Mo = 0.42 [IJ]) and a critical temperature T/J = 
1.89 [l5|. On the other hand, for g >> 1 the strong bonds 
form spin-singlet dimers, leading to the VBS phase. A 



QPT separates these phases, as shown schematically in 
Figfflb). This model has been studied previously |16j in 
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FIG. 4: (a) The model, with thin lines denoting J bonds and 
thick lines denoting gJ bonds; (b) Schematic phase diagram 
of the model at T = 0. 

connection with magnetic-field induced QPT's, and the 
quantum critical point was located at g^ = 4.013 ±0.003, 
using quantum Monte Carlo (QMC) methods. However, 
the important questions of the universal behaviour of 
the Neel temperature and the dynamics of the dimerized 
phase were not discussed. 

Our numerical calculations are based on series expan- 
sion methods [l7| . and involve several separate parts. 
The various series have been analysed in the usual way, 
via Pade approximants. The error bars shown on some of 
the data points are not statistical errors but 'confidence 
limits' based on consistency and spread between different 
high order approximants. For many data points, these 
error bars are smaller than the point size. 

We have used a 'dimer expansion' [l7[ to obtain series 
for the ground-state energy and the magnon energies in 
the VBS phase, in powers of l/p, to orders 11 and 8 
respectively. The latter provides a direct series for the 
minimum gap at k = (tt, tt, 0). Analysis of the gap series 
has to allow for the expected square-root singularity at 
gc, and we have used a Huse transformation to remove 
this singularity. The resulting gap data are shown in 
Fig. [1] by blue triangles. Our estimate of the critical 
point obtained from this data is fully consistent with, 
although somewhat less precise than the Monte Carlo 
estimate g^ = 4.013. We use this value in our further 
analysis. The gap data can be very well fitted by the 
expression 0.316^^(7 — pc, which is shown in Fig. [T]by the 
blue solid line. This provides the estimate A = 0.316Jg, 
where 

Jg = J(l+.9) 

is an average exchange parameter, used hereafter to set 
an energy scale. 

Results for the magnon energies near k = (7r,7r, 0), 
fitted to the expression 



efk) 



A2 + C2(^ - fci)2 + c2(7r - k^f -f c2fc| 



provide estimates of the magnon velocities near the QCP, 
c\ = C2 = 0.516Jg,c;3 = 0.337Jc,. These values contain 
an uncertainty up to ±5%. 

Next, we have used 'Ising expansions' [13] in the Neel 
phase to obtain scries for the ground state energy and 
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magnetization to order 12 in an anisotropy parameter 
X. These energy series, evaluated at x via Pade approx- 
imants, provide the data shown in Fig. 5 by red squares. 
The energies in the VBS phase, discussed above, are 
shown as black circles. As can be seen, the two en- 
ergy curves, from the Neel and VBS phases respectively, 
meet smoothly at the QCP, as expected for a second- 
order transition. This data can then be used to esti- 

-0.31| . ^ . ^ . ^ . 1 

-0.32- 




FIG. 5: The ground-state energy in the Neel phase (red 
squares) and dimer phase (black circles) versus the coupling 
constant. The lines are fits to the energy, as discussed below. 

mate the parameter a in Eq.(l4|). The VBS energy can 
be accurately fitted with a straight line EvBs/Jg = 
— 0.3244 — 0.0040(g— gc), shown as the black dashed line 
in Fig. [SJ The Nccl data can be fitted with a quadratic 
expression, as in Eq.dH), E^/Jg = EvBs{g) + 0.00003 - 
0.0101(g — QcY ■ However this fitting is subject to uncer- 
tainty, as the energies near the QCP are only changing 
in the 4th figure, and the data are not that precise. In- 
deed, inclusion of a small cubic term in the fit changes the 
coefficient of the quadratic term significantly, Ej^/Jg = 
EvBsig) - 0.00028 - 0.0088(5 - 9cY + 0.0008(<? - g^f. 
Comparing the coefficient of the quadratic term with 
Eq.ljl]) we determine the value of the quartic coupling 
constant. Our final estimate is a = 0.283Jg, with an 
uncertainty of ±15%. 

Substitution of the determined parameters into Eq.([7]) 
gives Tn/ Jg = 0.275-\/(?c — <?■ This dependence is shown 
in Fig. [T]by the solid magenta line with error bars. The 
main uncertainty, 7%, in the coefficient 0.275 comes 
from the uncertainly in the value of a discussed above. 
An additional few per cent come from uncertainties in 
A and magnon velocities. Altogether we estimate the 
computational uncertainly in the value of the coefficient 
0.275 as 10%. This is shown as error bars in the solid 
magenta curve in Fig. [TJ 

Finally, we compute high-temperature expansions for 
the Neel susceptibility. This is the response to a 'stag- 
gered' field. This susceptibility is expected to have a 
strong divergence at the critical temperature and can be 
used to estimate Tjv [g) ■ The Neel temperature calculated 
in this way is shown in Fig. [T] by the red squares. The 



red dashed line just connects the data points for guid- 
ance. The agreement between prediction of the universal 
theory shown by the magenta curve and results of the 
series computations is quite satisfactory. 

In summary, we have shown that a 3-dimensional an- 
tiferromagnet in the vicinity of a quantum critical point 
is expected to show universal behaviour, including scal- 
ing of the Neel temperature with the ground state mag- 
netization and with the coupling constant. We predict 
the universal scaling. Our prediction based on a field 
theory accurately describes recent data on the material 
TlCuCl^. The universal prediction is supported by nu- 
merical results obtained for a microscopic S=l/2 Heisen- 
berg spin model with strong and weak bonds, which is 
a specific example of a 3D antiferromagnet with a QCP. 
Results are obtained via a variety of series-expansion cal- 
culations, and are shown to be in reasonable agreement 
with the predicted universal behaviour, within numerical 
uncertainties. 
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